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J? ' Abstract. We present evolutionary synthesis models for galaxies of spectral types Sa through Sd with starbursts 

^^ , of various strengths triggered at various redshifts and study their photometric evolution before, during, and after 

^^ ' their bursts in a cosmological context. We find that bursts at high redshift, even very strong ones, only cause a 

\^ ' small blueing of their intrinsically blue young parent galaxies. At lower redshift, in contrast, even small bursts 

cause a significant blueing of their intrinsically redder galaxies. While the burst phase is generally short, typically 
a few hundred Myr in normal-mass galaxies, the postburst stage with its red colors and, in particular the very 
^ I red ones for early bursts at high redshift, lasts much longer, of the order of several Gyr. We find that, even 

^^ . without any dust, which in the postburst stage is not expected to play an important role anyway, models easily 

reach the colors of EROs in the redshift range z ~ 2 through z ~ 0.5 after starbursts at redshifts between 2 
and 4. We therefore propose a third alternative for the ERG galaxies beyond the two established ones of passive 
galaxies vs. dusty starbursts: the dust-free post-(strong-)starbursts. A very first comparison of our models to 
\^^ ' HDF data with photometric redshifts shows that almost all of the outliers that could not be described with our 

f^ , chemically consistent models for undisturbed normal galaxy types E through Sd can now be explained very well. 

Galaxies in the redshift range from z~ 2.5 to z~ 0.5 that are redder, and in some cases much redder, than our 
reddest undisturbed model for a high-metallicity classical elliptical are well described by post-starburst models 
after starbursts at redshifts between 2 and 4. Galaxies bluer than our bluest low metallicity Sd model, most of 
Q ' which have redshifts lower than 1, are well explained by ongoing starbursts. 

. . \ Key words. Galaxies: formation & evolution, starbursts, the Antennae (Whitmore & Schweizer 1995), with their 

P*" ' high-redshift global starbursts, or, even stronger, the Ultra Luminous 

k> IR Galaxies (ULIRGs), all of which feature extremely 

5^ - I , . strong nuclear starbursts, sometimes around an AGN, and 

5t . ' are in advanced stages of merging (cf. Genzel et al. 2001). 

1.1. Importance of starbursts at high redshifts a/t- -j-i, <-• ^ o i j +i, 

'^ ° Mmor mergers with mass ratios ^o.l under the same 

Despite considerable effort and advances from both the condition can trigger smaller starbursts. The situation is 

theoretical and observational sides, some very fundamen- less clear for accretion. Slow and steady accretion need not 

tal questions about the formation and evolution of galax- trigger a starburst directlybut may do so in an indirect 

ies are still open. One of them concerns the respective way by slowly destabilizing a disk. Accretion of wholesale 

roles of hierarchical merging versus accretion in the evo- lumps of gas (and stars) up to dwarf galaxies, on the other 

lutionary history of galaxies and their respective impacts hand, may enhance SF significantly. The boundary, how- 

on the star formation (SF) histories and on the final - ever, between enhanced SF and what would be called a 

or better today's - morphological and spectral appear- burst is not very precisely defined. Classically, a starburst 

ance of galaxies. Major mergers, i.e. mergers involving is defined as a rate of SF that is too high to be sustained 

two objects of comparable mass, clearly trigger strong for a long time without exhausting a galaxy's gas reser- 

starbursts, as long as they involve sufficient quantities of voir. With infall, this definition gets somewhat weakened, 

gas, as shown by examples in the local universe, like NGC ^^^^ ^i^j^i^^ ^^^ framework of a given structure and 

7252 (Schweizer 1982, Fritze - v. Alvensleben & Gerhard g^j^^^ formation scenario, fikc the hierarchical merging 

1994a fc b) or the Lum inous IR Galaxies (=LIRGs), hke p-^^^^.^^ ^^ ^^^^ extreme, and the initial collapse plus suc- 

Send offprint requests to: U. Fritze - v. A., e-mail: cessive accretion scenario - e.g. from filaments, as pro- 

ufritze@star.herts.ac.uk posed by Bournaud et al. (2005) - On the other, it is not 
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clear how much happens at which redshift. Recent de- 
tections of unexpectedly massive, evolved, and numerous 
galaxies at high redshift (Chapman et al. 2004, Cimatti 
et al. 2004, Franx et al. 2003, Le Floc'h et al. 2004) have 
challenged at least the semi-analytic version of the other- 
wise widely accepted A— CDM models. 

Many different and independent pieces of evidence in- 
dicate that the merger rate was higher in the past. How 
much higher, however, is still controversial (e.g. Zepf & 
Koo 89, Wu & Keel 1998, Consehce 2003). 

Galaxies were richer in gas in the past, so that mergers 
should have caused even stronger starbursts than locally. 
Damped Lyman-a systems, abundant in number density 
at z ^ 2 — 4 and containing the bulk of the baryonic mat- 
ter at these redshifts in the form of HI (e.g. Wolfe et al. 
1993), are beheved to be progenitors of present-day spi- 
rals and to already contain tremendous amounts of mass, 
almost all of it still in the form of gas (Prochaska & Wolfe 
1998). If those objects were to be involved in a merger, 
the reservoir for a starburst would be gigantic. Indeed, 
there is even direct evidence that starbursts are stronger 
towards higher redshifts. The space densities of LIRGs 
and ULIRGs increase with redshift, and new classes of 
extremely powerful starbursts appear: hyper-luminous IR 
galaxies and SCUBA sources. Their dust-enshrouded SF 
rates are estimated to be of the order of thousand(s) of so- 
lar masses per year (e.g. Chapman et al. 2004, Le Floc'h 
et al. 2004). 

Many Lyman Break Galaxies (LBGs) are also called 
starbursts by some authors, albeit with more moderate 
SF rates than ULIRGs or SCUBA galaxies, typically in 
the range 10 - 100 Mg/yr (Shapley et al. 2001, Pettini 
et al. 2001, Giavalisco 2002). Although some of them fea- 
ture irregular or knotty structures, they are not clearly 
related to (major) mergers. It should be kept in mind 
that even the good old monolithic initial collapse sce- 
nario (Eggen, Lynden - Bell & Sandage 1962), as well as 
many different pieces of evidence - from the present stel- 
lar masses through chemical enrichment to fossil records 
- predict, that the SF rates of normal Hubble sequence 
galaxies should have been higher in the past (e.g. Sandage 
1986). Hence, the discrimination between "normal SF" 
and a burst also gets less clear-cut towards larger look- 
back times. 

All in all, the currently emerging scenario of global 
cosmic SF describes two phases: at high redshift SF in 
galaxies was dominated by discrete, recurrent bursts, pos- 
sibly associated with mergers. Around z ^^ 1.4 the universe 
switched gears and SF in galaxies became more quiescent, 
governed by accretion more than by mergers, allowing the 
Hubble sequence to emerge (cf. Papovich et al. 2005). 
With respect to individual galaxies, both the mass assem- 
bly and the bulk of SF appear to have occurred earlier in 
galaxies that by today are very massive than in lower mass 
galaxies where the bulk of present-day SF is still taking 
place. Extending results by Heavens et al. (2004) on the 
fossil record of SF in ^ 10^ nearby galaxies from the SDSS, 
Hammer et al. (2005) present additional evidence from a 



sample of CFRS galaxies at z > 0.4 for a scenario in which 
intermediate mass galaxies (3 — 30) x 10"'^° M© formed a 
large fraction of their stars only 4 to 8 Gyr ago, most prob- 
ably in a series of recurrent bursts of SF, many of which 
apparently gave rise to LIRG phases. The high fraction of 
LIRGs among intermediate mass galaxies at z > 0.4, to- 
gether with a fraction of > 17% of those LIRGs showing 
evidence for major mergers, supports this hypothesis. 

All this clearly shows the important role of starbursts 
(and mergers) at high redshift. 

1.2. Motivation for the study of starburst and 
post-starburst models 

A comparison between our set of undisturbed, chemically 
consistent galaxy models E through Sd over the redshift 
range from z > 4 through z > and the Hubble Deep Field 
North (HDF-N) galaxies with photometric redshifts from 
Sawicki et al. (1997) was presented in Bicker et al. (2004). 
We recall that the agreement between our chemically con- 
sistent models, which account for the increasing metal- 
licities of successive stellar generations, and the data was 
much better than for any earlier models using solar metal- 
licity input physics alone. The inclusion of stellar subpop- 
ulations with lower-than-solar metallicity made the chem- 
ically consistent Sd model bluer at increasing redshifts 
than an equivalent model with solar metallicity alone. On 
the other hand, including stars with supersolar metallicity 
made our chemically consistent, classical elliptical model 
redder at intermediate redshifts than the corresponding 
model for Zq. Nevertheless, while we found our chemically 
consistent galaxy models to bracket the redshift evolution 
of the bulk of the HDF data well over the redshift range 
from z > 4 through z > between the red elliptical and 
the blue Sd models, there remained a number of galax- 
ies with bluer colors than those of our bluest undisturbed 
Sd model, described by a constant SFR and fairly low 
metallicity, at low and intermediate redshifts, as well as 
quite a few galaxies with significantly redder colors than 
for our reddest undisturbed elliptical galaxy model - de- 
spite its high metallicity in the redshift range 0.5 < z < 4. 
This motivated us to investigate if and to what extent 
starbursts can explain the "bluer-than-normal" galaxies 
and if and to what extent post-starbursts can explain the 
"redder-than-normal" ones. 



1.3. The present study 

We here present models for galaxies with starbursts of var- 
ious strengths occurring at various evolutionary stages or 
redshifts to study the impact of ongoing and past star- 
bursts on the photometric properties of these galaxies. 
While spectra, of course, can give fairly precise SF rates, 
the advantage of the photometric signatures of starbursts 
and post-starbursts is that they can be compared to the 
wealth of multi- wavelength imaging data available for sev- 
eral deep fields, and help constrain the rates, strengths. 
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and fate of starbursts for huge numbers of galaxies, most 
of which arc too faint for spectroscopy. 

The widely used SED (spectral energy distribution) fit- 
ting technique to obtain photometric redshifts and galaxy 
types - using observed or model SEDs - gives the num- 
ber of ongoing starbursts as a function of redshift. Only 
the evolutionary synthesis approach presented here, how- 
ever, allows a galaxy to be followed through a starburst 
and beyond in a cosmological context and, hence, to iden- 
tify the successors of high-redshift starbursts in terms of 
post-starburst galaxies at lower redshifts. 

The outline of the present paper is the following. 
In Sect. 2.1, we briefly present our evolutionary synthe- 
sis code GALEV that has already been extensively de- 
scribed in the context of evolution models for star clusters 
and "normal" undisturbed galaxies, including cosmolog- 
ical and evolutionary corrections from high to low red- 
shift. We describe how we put starbursts on top of our 
undisturbed galaxy models and discuss the grid of pa- 
rameters we chose to study in Sect. 2.2 and 2.3. In Sect. 
3.1, we present selected results and show the impact of 
bursts of various strength at various times/redshifts on 
the luminosities and colors of galaxies from the optical to 
NIR, not only during the active burst phase but also af- 
ter the bursts. In Sect. 3.2 and 3.3, we discuss different 
burst strengths and the role of dust. In Sect. 3.4, we in- 
vestigate in how far our models relate to the Extremely 
Red Objects (EROs) in their postburst phases and pro- 
pose a new scenarion for the passively evolving 50 % of 
the ERO population in terms of (dust-free) post-strong- 
starbursts. Finally, Sect. 3.5 shows a first comparison with 
Hubble Deep Field (HDF) data. A consistent interpreta- 
tion of HDF galaxies in terms of our chemically consistent 
models for undisturbed galaxies, the starburst and post- 
burst models as presented here, our own photometric red- 
shifts and a detailed quantitative analysis in terms of burst 
strengths, starburst and post-starburst galaxy fractions, 
metallicity and mass evolution, etc., will be the subject of 
a forthcoming paper. 



2. Models 

2.1. Undisturbed galaxies 

Our evolutionary synthesis code GALEV in its latest ver- 
sion has been described in the context of the chemically 
consistent evolution of galaxies (Bicker et al. 2004). The 
code accounts for the increasing initial metallicitics of suc- 
cessive stellar generations and naturally reproduces the 
observed broad stellar metallicity distributions of local 
galaxies. It is based on isochrones for 5 different metal- 
licities —1.7 < [Fe/H] < +0.4 from the Padova group that 
include the thermal-pulsing AGB phase, which has been 
shown by Schulz et al. (2002) to have an important ef- 
fect on colors like (V — I) and (V — K) at ages from 100 
Myr through a few Gyr for single-burst, single-metallicity 
stellar populations like star clusters. 



The spectral evolution of a stellar population is cal- 
culated on the basis of the Lejeune et al. (1997, 1998) 
library of model atmosphere spectra for the five metallic- 
ities. Luminosities in a wide variety of filter systems and 
colors are obtained by folding the spectra with the re- 
spective filter response functions. In Anders & Fritze ~ v. 
Alvensleben (2002) we added the emission contributions 
from gas ionized by young stars both in terms of contin- 
uous and line emission for the respective metallicitics to 
our single burst single metallicity models, while in Bicker 
& Fritze - v. Alvensleben (2005) we discussed its effects 
in the context of galaxies. 

The various spectral types of undisturbed galaxies 
from E through Sd are described by their respective ap- 
propriate SF histories. The spiral models Sa, Sb, Sc, Sd 
discussed in the present context use SF rates linearly tied 
to the evolution of the gas content that decreases through 
SF and increases through stellar mass loss. For the model 
of a classical initial-collapse elliptical galaxy, we use a SFR 
exponentially declining with time with an e-folding time of 
1 Gyr. Gonstants are chosen so as to match, after ^13 Gyr, 
the observed average colors, gas content, characteristic HII 
region metallicities, M/L-ratios, and template spectra of 
the respective spectral types in the local universe. Models 
do not include any spatial resolution, dynamics, or AGN 
contributions. No infall is included, galaxies are described 
by closed box models for simplicity. We have shown that 
continuous infall that increases the mass of a galaxy by no 
more than a factor of ^ 2 does not significantly affect the 
time evolution of spectra, colors, and metallicities, since 
SF histories in this case need to be adjusted in order to 
still match the present-day observables. To a large extent, 
this cancels the effects of this type of infall on the color 
and metallicity evolution (Fritze - v. Alvensleben 2000). 

A standard cosmological model with 

Hq — 70, ilm — 0.3, i^A = 0.7, and a redshift of galaxy 
formation of Zf ^ 5 are assumed to calculate redshifted 
galaxy spectra, luminosities, and colors. For galaxies at 
high redshifts the average attenuation by intergalactic HI 
along the line of sight as determined by Madau (1995) 
is included in the models (cf. MoUer et al. 2001). Dust 
within the galaxies, however, is not yet included at the 
present stage, neither in terms of absorption in the UV 
and optical nor in terms of thermal recmission in the 
mid- and far-IR. 

Gaseous emission in terms of lines and continuum 
is consistently included in our models for the respective 
metallicities (cf. Bicker & Fritze - v. Alvensleben 2005). It 
has been shown to have a large impact, even on broad 
band luminosities and colors, during strong starbursts 
by Kriiger et al. (1995) (see also Anders & Fritze - v. 
Alvensleben 2003). 

Our undisturbed model galaxies are scaled to match 
the mean present day observed Mb of the respective 
galaxy types in the local universe as determined by 
Sandage et al. (1985). These average absolute B-band lu- 
minosities and their 1-cr ranges are (Mb) = — 20.8 ± 1.7, 
(Mb) = -18.1 ± 0.94, and (Mb) == -17.0 ± 0.77 mag for 
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Fig. 1. Example of the SFH for a bursting galaxy. In this 
case a strong burst starts in an Sb galaxy at an age of 4.7 
Gyr, corresponding to a redshift z = 1. Top panel: time 
evolution, bottom panel : redshift evolution. 

E, Sb, and Sd galaxies, respectively (cf. Bicker et al. 2003 
& 2004). 

2.2. Starbursts 

As a consequence of the assumption of ideal and instan- 
taneous gas mixing and for the closed-box approximation 
used here, we cannot treat starbursts in a chemically con- 
sistent way. We therefore chose to model the bursting 
galaxies with fixed half-solar metallicity input physics in 
fair agreement with the high redshift galaxies with spec- 
troscopic abundance determinations (e.g. Mehlert et al. 
2002, Shapley et al. 2005). For consistency, the undis- 
turbed galaxies used for comparison in this paper are also 
calculated with fixed half solar metallicity. 

The bursts are described by a sudden increase in the 
SFR at the onset of a burst at redshift Zburst or time tburst 
to a value ^max, followed by an exponential decline with 
an e-folding time or burst duration Tburst (cf- Fig. 1). 

This simple scenario is seen to cope well with obser- 
vations of local starbursts and post-bursts that become 
passive after the end of the burst due to gas consump- 
tion by SF and, eventually, a loss of left-over gas in a 
starburst-driven wind. We caution, however, that in the 
case of starbursts at high redshifts a galaxy may regain 



some SF activity after a while, either due to accretion of 
primordial gas or due to re-accretion of gas that was for- 
merly expelled. 

The strength of the burst b is defined by the fraction 
of gas available at the onset of the burst G (tburst), which 
is consumed for star formation during the burst (AG): 

Another common definition of the burst strength uses the 
fraction of stars formed in the burst as compared to the 
total stellar mass of the galaxy. We caution that, due to 
the evolving stellar-to-gaseous mass ratio, there is no 1- 
to-I relation between these two definitions. 

For clarity, we separately study starbursts within a 
given galaxy and starbursts induced by major (= 1:1) 
mergers. All our figures refer to the first class of star- 
bursts within a given galaxy. A major merger between 
two galaxies of comparable type and luminosity leads to 
an additional increase in luminosity by 0.75 mag without 
affecting the colors. 

2.3. Parameter ranges to be explored 

Our undisturbed spiral galaxy models Sa, Sb, Sc, Sd are 
normalized to give by today, i.e at an evolutionary age of 
to ~ 13 Gyr, the respective mean local B-band absolute 
luminosities observed by Sandage et al. (1989) for the re- 
spective galaxy types. Bursts are assumed to be triggered 
at a range of redshifts z = 4, 3, 2, 1, 0.5. Burst durations 
were shown by Fritze - v. Alvensleben & Gerhard (1994a) 
to be a parameter that cannot be very tightly constrained. 
It is to some degree degenerate with the detailed shape of 
the burst, i.e. steeply/slowly rising/falling SFRs. Hence, 
we chose a fixed burst duration, i.e. e-folding time for the 
exponential decrease in the SFR of t;, — 2.5 x 10^ yr, as 
defined in Sect. 2.2. In general, burst durations appear to 
be related to the dynamical timescales in galaxies. While 
dwarf galaxies feature burst durations on the order of 10^ 
yr (Kriiger et al. 1991), normal galaxies - the objects of 
the present study - are generally described by burst dura- 
tions on the order of a few 10® yr (e.g. Barnes & Hernquist 
1996, Fritze - v. Alvensleben & Gerhard 1994a). 

The burst strengths we explore for this study range 
from a consumption of 20% to 50% and 80% of the avail- 
able gas reservoir at the onset of the burst for weak, inter- 
mediate, and strong bursts, respectively. A general one- 
to-one transformation into the relative increases in stellar 
mass is not possible since the stellar mass and the mass of 
the gas reservoir vary with redshift and galaxy type. For 
example, an intermediate strength (50%) burst on top of 
an Sb type galaxy will increase the stellar mass by a factor 
of 8.6 for a burst at redshift z = 4, by a factor 2.3 for a 
burst at z = 2, and by a factor 1.2 for a burst at z = 0.5. 

Note that at the present stage, our models do not in- 
clude any dust, wether in absorption or in its thermal 
reemission at rest-frame IR wavelengths. A consistent in- 
clusion of appropriate amounts of dust in the course of the 
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redshift evolution of the various model galaxy types is dif- 
ficult, in particular during bursty phases of SF. Optical as 
well as FUV detected starbursts in the local universe along 
with LBGs at high redshift, all do show dust extinction, al- 
beit typically at rather moderate values of Apuv '^0.5 — 2 
mag (Buat et al. 2005, Shapley et al. 2001, Pettini et al. 
2001), corresponding to visual extinctions on the order 
of Ay ^ 0.2 — 0.9 mag. A comparison between FIR (60 
/im) and NUV (GALEX) selected starburst galaxy sam- 
ples shows that all but one out of 118 FIR-selected star- 
bursts are also detected in the NUV and that the median 
NUV dust attenuation for the FIR-selected starbursts is 
only ^ 2 mag (Buat et al. 2005). Hence, while the present 
models clearly cannot be applied to the interpretation of 
very dusty starbursts like those observed in ULIRGs or 
SCUBA galaxies, they will certainly give a valid approxi- 
mation of the majority of more normal, i.e. not too dusty, 
starbursts. 

3. Results 

3.1. Impact of starbursts on the photometric 
properties of galaxies 

For our discussion of the model properties we use photo- 
metric data in the standard Johnson-Cousins filters in the 
VEGAMAG system. 

In Fig. 121 we show the B-band luminosity evolution as 
a function of galaxy age on the example of an Sb galaxy 
with intermediate-strength bursts at various times corre- 
sponding to redshifts 4 > z > 0.5. After an initially rapid 
increase, the luminosity of the undisturbed Sb galaxy 
reaches a maximum around 2 Gyr. In the following 10 
Gyrs its absolute luminosity Mb fades by ~ 1 mag. During 
a burst the galaxy gets significantly brighter, as expected. 
Due to the exponential decline of the burst SFR on a 
timescale of 2.5 x 10* yr, the luminosities after the bursts 
fade rapidly below the level of the undisturbed Sb galaxy. 
The bursts in Fig.|21all have the same strength of b = 50% 
in our definition of gas consumption. At younger ages, 
when the galaxy has more gas available, the luminosity 
increase is higher for a given burst strength than at older 
ages, when the gas reservoir is smaller. 

We clearly see that none of the bursts gets as bright 
as a classical initial collapse elliptical model (the lowest 
curve in Fig. 2, upper panel) plotted for comparison. Long 
enough {'^ 1/2 tnubbic) after early bursts, the Sb post- 
burst model gets fainter, without any further SF, in Mb 
than even the faintest and latest-type undisturbed galaxy 
model Sd. We recall that our undisturbed model galax- 
ies are scaled to match the mean present day absolute 
B-band luminosities (Mb) of the respective galaxy types 
in the local universe. 

The time evolution of the (B — V) color shows that 
bursts occurring at later phases, although transforming 
smaller absolute amounts of gas, cause stronger blueing, 
since by that time the galaxy itself is already relatively red 
than bursts occurring at earlier stages, when the galaxy 
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Fig. 2. Time evolution of the absolute B-band lumi- 
nosity Mb (top panel) and of the restframe (B — V) 
color (bottom panel) of an Sb galaxy with intermedi- 
ate strength bursts at times corresponding to redshifts 
z = 4, 3, 2, 1, and 0.5. For comparison, the undisturbed 
Sd and E models are also displayed to mark the upper 
and the lower luminosity limits and the red and blue color 
extremes of our undisturbed galaxy models. 



itself is still young and blue. On the other hand, the red- 
dening after the burst is much stronger during the long 
evolution after early bursts than during the shorter time 
interval between a more recent burst and the present. A 
color of (B — V) ~ 0.9 is a red limit to any old, passively 
evolving stellar population. Within a Hubble time galax- 
ies cannot get any redder than that, even after > 6 Gyr 
of passive evolution without any SF. 

It is seen that during and up to ^ 600 Myr after the 
bursts, the Sb-|-burst models get significantly bluer than 
an undisturbed Sd. Shortly thereafter, however, the post- 
burst models reach colors as red ~ or even much redder 
after early bursts - than those of a passively evolving el- 
liptical. Note that while the blueing during a burst is a rel- 
atively short-time phenomenon lasting < 300 — 500 Myr, 
the red and very red phases after the end of a burst can 
last very long, i.e. for many Gyr, unless a galaxy resumes 
SF from gas (re-)accreted after the end of the burst. 

Figure 3 shows the evolution of the same models as a 
function of redshift. The luminosity evolution, however, is 
now given in terms of observer-frame apparent Johnson 
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Fig. 3. Redshift evolution of an Sb galaxy with intermedi- 
ate strength bursts at z = 4, 3, 2, 1, and 0.5 in terms of 
observer-frame, apparent B-band luminosities me (upper 
panel) and in terms of observer- frame, apparent Ic— band 
luminosities mi^ (lower panel). Undisturbed Sd and E 
models are again displayed for comparison. 



B-band and Cousins I-band luminosities (upper and lower 
panels, respectively). Now, the luminosity evolution is no 
longer only an evolutionary effect, but also includes the 
cosmological bandshifts, dimming, and attenuation. The 
behavior of the bursting galaxies is similar to that ex- 
plained above. Note, however, the very strong apparent 
fading in niB of galaxies after starbursts at high redshift 
by the combined effects of fading and bandshift; e.g., an 
Sb galaxy with an intermediate strength starburst at z ~ 3 
gets brighter by about 2.5 mag in niB until a redshift 
z ~ 2.5 and then becomes very faint ~ almost 4 mag fainter 
in niB than it was before the burst and almost 5 mag 
fainter by a redshift around 1.5 than an undisturbed Sb 
galaxy. The same galaxy with the same starburst around 
z ~ 2 would be about 4 mag fainter by z ^ 1 than an 
undisturbed galaxy. 

Another effect results from the nonlinear transforma- 
tion between time and redshift. While the e- folding time 
is the same for all bursts (2.5 x 10^ yr), the duration - 
in terms of redshift - varies. A burst occurring at z = 4 
will light up the galaxy until z = 3. In contrast, a burst at 
z = 1 will light up the galaxy only up to a redshift z = 0.8. 



Fig. 4. Redshift evolution of Sb models with intermediate- 
strength starbursts at different redshifts in observer-frame 
B — V (upper panel) and V — Ic (lower panel) . 

As seen in the lower panel of Fig. 3, the effects of 
starbursts on apparent I-band luminosity are qualitatively 
similar to those in B, but much smaller both in terms of 
brightening during and fading after the burst. 

Figure 4 shows the redshift evolution of (B — V) and 
(V — Ic) colors in observer- frame. This is more com- 
plicated to understand than the time evolution in Fig. 
2. Both bandshift effects and attenuation play a role. 
Those Sb galaxies with bursts at redshifts > 3 get bluer 
not only during the active burst but also in the post- 
burst phase until redshifts z = 2 — 2.5, in line with the 
undisturbed Sb models. Nevertheless, these postbursts 
are always slightly redder than the undisturbed galaxy. 
Bursts at z = 4, z = 3, and z = 2 become very red, up 
to (B — V) '^ 2 by redshifts z = 2 , z = 1.5, and z = 1, re- 
spectively, through z ~ 0.5. In terms of (V — Ic), these 
same postburst galaxies can reach values as red as 
(V — Ic) = 2 — 3.2 over the redshift range 2 > z > 0.5. 

In these color vs. redshift plots (Fig. 4), it is clearly 
seen that the reddening after starbursts at high redshifts 
(4 > z > 0.5) has a very strong effect, one much stronger 
than the blueing during the burst. Comparison with the 
classical initial collapse model for ellipticals also plotted in 
these figures shows that, provided SF stops completely, the 
post-starburst models reach colors at redshifts between 
0-5 "^ z <2 that are redder than those of the passively 
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evolving E-model by 0.5 to 2.0 mag in (B — V) and by 1 
to 1.5 in (V — I). The postburst models reach color ranges 
never attained by undisturbed or passively evolving galax- 
ies: (B - V) > 1.5 and (V - Ic) > 1.5 at 0.5 < z < 2 for 
bursts at 2 < z < 4. We caution that these red and faint 
postburst stages will only be reached by galaxies that do 
not resume any residual SF after a burst, e.g. from previ- 
ously expelled and later re-accreted enriched gas or from 
newly accreted unenriched material. The earlier a burst 
occurs, the more time a galaxy has until today to reassem- 
ble a gas reservoir for subsequent SF, as discussed in the 
context of secondary disk rebuilding (e.g. Hammer et al. 
2005). On the other hand, the time between a starburst at 
z ~ 4 (or 3) and the redshift of maximum postburst red- 
dening at z ~ 2 ... 0.5 (or z ~ 1.5 ... 0.5) is short. In the 
case of a burst at z ~ 4, the phase of maximum reddening 
starts 1.7 Gyr after the burst and holds on for 5.3 Gyr. 
After intermediate-strength or strong bursts, in particu- 
lar, which most probably lead to a strong galactic wind, as 
observed in many, if not all, high-redshift starburst galax- 
ies, the galaxies will most probably need more time than 
this before they will eventually be able to restart subse- 
quent SF. We thus expect to really observe passive post- 
starburst galaxies with very red colors (B — V) ~ 1 — 2 
and (V — Ic) ~ 1.5 — 3 in the redshift range of z ^ 2 to 
z - 0.5. 

Figure |S1 briefly highlights the spectral evolution of 
our models, from which the photometric properties of red- 
shifted galaxies are derived. The separate panels show a) 
the actively SFing Sb galaxy before the onset of the burst, 
b) the enhanced short-wavelength flux and the strong 
emission lines during the burst, c) and d), respectively, 
the model spectra 0.5 and 1 Gyr after the onset of the 
burst. Panel c) gives a glimpse of the strong higher order 
Balmer lines during the post-burst phase and d) shows 
the almost completed transformation from a SFing into 
a passive K-star-type spectrum with weak emission lines 
from the starving star formation still on top. 

3.2. Ef Facts of the burst strength 

Before, all the bursts we studied had the same strength 
b — 0.5, i.e. transforming 50% of the galaxy's remain- 
ing gas reservoir into stars. Now we study the influ- 
ence of different burst strengths. Figure 6 shows the red- 
shift evolution during and after bursts with strengths 
b = 0.2, 0.5, and 0.8 occurring in an Sb galaxy at a red- 
shift z = 1. In terms of apparent B-band luminosity, the 
differences are clearly visible, in particular between a weak 
and an intermediate strength burst, and less so between 
the intermediate and strong bursts (top panel). 

In (B — V), as in all optical colors, the differences be- 
tween the different burst strengths are fairly small (lower 
panel). Even smaller effects on the colors are observed in 
our models for bursts occurring at higher redshifts when 
the blueing effect is intrinsically smaller because of the 
bluer colors of the underlying undisturbed galaxy models. 
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Fig. 5. Time evolution of the optical spectra for an Sb 
model with an intermediate strength starburst occurring 
at a redshift z = 1: (a) prior to the burst, (b) at maximum 
burst, (c) 0.5 Gyr after the burst, and (d) 1 Gyr after the 
burst. 




Fig. 6. The impact of different burst strengths on the red- 
shift evolution of mi, (top panel) and on the color evolution 
(B — V) (lower panel) for starbursts occurring at z = 1 in 
an Sb galaxy and consuming 20%, 50%, and 80% of the 
available gas, respectively. 
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3.3. The role of dust 

Although our models at their present stage do not include 
any dust, we briefly want to discuss its effects. If we were 
to include realistic amounts of dust in our undisturbed 
models, we would have to slightly adjust the SFRs of the 
respective spiral models so as to again reproduce the av- 
erage observed colors after a Hubble time (cf. MoUer et 
al. 2001). Hence, the inclusion of typical, i.e. moderate, 
amounts of dust would not affect the evolution of colors 
by much. During the bursts, dust could significantly re- 
duce the blueing and the luminosity, in particular in short 
wavelength bands, if present in sufficient quantities, as 
in the case of dusty starbursts like LIRGs or ULIRGs in 
the relatively local Universe or SCUBA galaxies in the 
more distant one. Shortly (< 500 Myr) after the peak of a 
starburst, however, local example galaxies appear to have 
largely or even completely cleared out their dust (cf. NGC 
7252, NGC 3921, . . .). As a result, we expect dust to no 
longer play any significant role during most of the post- 
starburst phase, no matter how dusty the starburst phase 
was. 

3.4. Relation to EROs 

The nature and origin of the Extremely Red Objects, 
called EROs, defined by having (R — K) > 5, which are 
found in ample number densities at redshifts z ~ 1 — 2, is 
a topic still under some discussion. Recent results suggest 
that the ERO galaxy population is composed to approxi- 
mately comparable parts of dusty starbursts and old pas- 
sively evolving stellar systems. The passive galaxies are 
thought to be the high-luminosity high-mass end of the 
elliptical galaxy progenitor population at 1 ^ z < 2 and 
may be related to the starburst galaxy population at z > 2 
seen in deep submillimeter surveys (cf. Mc Carthy 2004 for 
a recent review). Moustakas et al. (2004) report ^ 40% of 
all EROs in the GOODS field to show early-type mor- 
phologies. 

Comparison with our models in terms of (R — K) (Fig. 
7, upper panel) indicates that the classical initial collapse 
model for an elliptical galaxy does not reach the color 
range of EROs. This is only valid, however, for an average 
luminosity elliptical. The luminosity-metallicity relation 
observed for elliptical galaxies leads us to expect that an 
elliptical galaxy that by today will be brighter than aver- 
age will also have higher metallicity than our model and, 
hence, reach redder colors and possibly go though an ERO 
stage in the course of its evolution - provided it formed in 
this classical scenario and evolved in isolation, i.e. without 
accretion of gaseous matter triggering star formation. Our 
postburst models with 2 < Zburst ^ 5, however, all reach 
far into the ERO color range between z r^ 2.3 and z ^^ 1. 

Comparison with the classical E-model shows that 
these very red colors are only achieved by stellar popu- 
lations that formed the bulk of their stars within a short 
time interval on the order of a few 10^ yr and rapidly 
dropped their SFRs to zero thereafter. The tail of the 1 




Fig. 7. Redshift evolution of Sb models with intermediate 
strength starbursts at different redshifts in observer frame 
R — K (top panel) and in mx. 



Gyr exponentially declining SFR in our classical E-model 
is sufficient to keep colors out of the extremely red (R — K) 
range forever. 

Another property of observed EROs is their high lu- 
minosity K > 20 mag. In the lower panel of Fig. 7 we plot 
the redshift evolution of apparent K-band magnitudes niK 
for our Sb-f intermediate strength burst models, as well as 
for our undisturbed E, Sa, and Sd models. The K-band 
(and other NIR) luminosities increase during bursts, al- 
though somewhat less than e.g. the B-band luminosity. 
An intriguing feature of NIR luminosities, in contrast to 
blue and optical ones, is that they do not significantly de- 
crease after the burst, but rather remain fairly constant 
instead. This is due to the well-known fact, that NIR light 
traces mass and that mass is dominated by long-lived low- 
mass stars, while blue and optical light instead traces the 
short-lived higher- mass stellar component. Looking at the 
K-band evolution of our models in the lower panel of Fig. 
7, we note that our intermediate strength bursts in an 
average-luminosity Sb galaxy are not luminous enough to 
reach the K< 20 domain. However, if a burst, and in par- 
ticular a strong one, takes place in a more luminous early 
type spiral or is triggered by a major merger that doubles 
the luminosity, the K< 20 criterium of EROs is easily met 
by galaxies with starbursts at z < 4. 
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As discussed in Sect. 13.31 dust should not be an is- 
sue during most of the post-starburst evolution, so we do 
expect our dust-free models to give approximately cor- 
rect colors and luminosities in this stage. Even if occur- 
ring in a spiral-type galaxy, we expect a reasonably strong 
burst - which, at high redshift, is most probably triggered 
by interactions and mergers during hierarchical structure 
formation - to be accompanied by a restructuring of the 
remnant towards an E/SO morphology type, in agreement 
with the observations by Moustakas et al. (2004) men- 
tioned above. While our dust-free models do not give any 
clue to the dusty-starburst part of the ERO galaxy pop- 
ulation, they show that the ~ 50 % so-called passively 
evolving stellar systems are explained well by our dust- 
free post-starburst models and not by average-luminosity 
classical E-models with an initial collapse and subsequent 
passive evolution. Our results indicate, that the K < 20 
objects require strong - probably merger-induced - star- 
burst progenitors. While brighter-than-average, classical, 
passively evolving elliptical models might work due to the 
color-luminosity relation, we propose here a new scenario 
for ERO galaxies, the dust-free post-strong-starbursts, as 
a third alternative to the conventionally assumed, pas- 
sively evolving vs. dusty-starburst dichotomy. 



3.5. First comparison witti HDF data 

The incentive for the present investigation of the role of 
starbursts on the photometric evolution of high redshift 
galaxies came from a comparison of our chemically consis- 
tent evolutionary synthesis models for undisturbed galax- 
ies with data for intermediate and high redshift galaxies 
in the Hubble Deep Field North (HDF-N) (cf. Bicker et al. 
2004). While we found good agreement between our chem- 
ically consistent galaxy models and the bulk of the HDF 
data over the redshift range from z > 4 through z > 0, 
there remained a number of galaxies at intermediate and 
high redshifts with bluer colors than in our bluest undis- 
turbed Sd model described by a constant SFR and fairly 
low metallicity, as well as quite a few galaxies with col- 
ors significantly redder in the redshift range 0.5 ;^ z < 4 
than our reddest undisturbed elliptical galaxy model, de- 
spite its high metallicity. This lead us to investigate to 
what extent starbursts can explain the bluer than normal 
galaxies and post-starbursts the redder than normal ones. 

Therefor, we now compare our starburst and post- 
starburst models with photometric data for galaxies in 
the HDF using the photometric redshifts given by Sawicki 
et al. (1997). For this comparison, we calculate the cor- 
responding photometry from our model spectra explicitly 
for the HDF filters in the AB-mag system to avoid uncer- 
tainties involved in a transformation from HST to Johnson 
filters. In Fig.|Slthc observed m450 (^ B) magnitudes are 
plotted together with the redshift evolution of our burst 
and post-starburst models and with our undisturbed E, 
Sb, and Sd models for comparison. The bars in Fig.|Hlshow 
the Icr— luminosity ranges for the respective undisturbed 
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Fig. 8. Redshift evolution of our models in terms of HST 
450W apparent luminosities mm45o compared to HDF-N 
galaxies with photometric redshifts from Sawicky et al. 
(1997). 




Fig. 9. Color evolution in terms of (m450— m606) for star- 
burst and undisturbed models compared to HDF-N data. 



local galaxy types. We caution that these may change with 
redshift due to bandshift effects. 

The undisturbed model for the classical elliptical 
galaxy is much brighter than any observed galaxy at red- 
shifts greater than z = 2. This indicates that the classical 
monolithic initial collapse model is not a likely scenario 
for today's elliptical galaxies. More than half of the HDF 
galaxies at 0.5 < z < 1.2 are fainter in m450 (^ B) than 
our faintest galaxy model Sd. The brighter of these could 
be actively SFing dwarfs, while the fainter ones in any 
case must be passive galaxies with very little emission in 
rest-frame U. At the lowest redshifts z < 0.5 the selection 
criteria for the HDF result in a lack of bright galaxies. 
Sb models with intermediate strength bursts easily reach 
the brightest galaxies in the HDF-N at all redshifts during 
their active burst phase, as well as almost all of the faint 
galaxies at redshifts between z = 1.5 and z = 0.5. 

More information is obtained from a comparison be- 
tween observations and model galaxies in terms of col- 
ors. In Fig. 1^1 the redshift evolution of (m450— m606) 
~ (B — V) is plotted. While the bulk of the observed 
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galaxies fall well within the color range of our undisturbed 
galaxy models, some are bluer than even our bluest Sd- 
model over the entire redshift range from z > 4 through 
z ~ and some are significantly redder than our undis- 
turbed E-model that marks the red envelope of our undis- 
turbed models. Around a redshift of z = 1 there are a 
number of galaxies with (m450 — m606) > 1.5 that can- 
not be reached by our undisturbed E model, but only 
by post-starburst galaxies after bursts at redshifts z > 2. 
Essentially all galaxies in the redshift range 1.5 < z < 3.5 
with colors much redder than those of a classical passively 
evolving elliptical model are reached by post-starburst 
models. In the redshift range 1 ^ z < 2, there seems to be 
an apparent lack of galaxies with the colors of our post- 
burst models. If the red galaxies around z = 0.5 ... 1 are 
supposed to be postburst galaxies, there should also be 
some postburst galaxies in this area. A look back at the 
luminosity evolution (cf . Fig. |SJ solves this apparent prob- 
lem: in this redshift range the postburst galaxies are too 
faint and fall below the detection limit for the HDF. All 
HDF galaxies bluer than our bluest undisturbed Sd model 
are easily reached by ongoing starbursts and at low red- 
shift, even a weak burst can cause a significant blueing as 
shown in Figs. 4 and 6. 

A detailed and quantitative analysis of the HDF 
galaxy population in terms of undisturbed, starburst, 
post-starburst models, SFRs, metallicities, masses, etc., 
using our own consistent photometric redshift determina- 
tions is currently under way and will be the subject of a 
forthcoming paper. 

4. Conclusions and outlook 

The incentive for this present investigation came from a 
comparison of our chemically consistent evolutionary syn- 
thesis models for undisturbed galaxies with data for in- 
termediate and high redshift galaxies in the Hubble Deep 
Field North (HDF-N) (cf. Bicker et al. 2004). The agree- 
ment of our chemically consistent models, which account 
for the increasing metallicities of successive stellar gener- 
ations, with the data was much better than for any earlier 
models solely using solar metallicity input physics. The 
inclusion of stellar subpopulations with lower-than-solar 
metallicities made the chemically consistent Sd model 
bluer at increasing redshifts than an equivalent model with 
only solar metallicity. The inclusion of stars with super- 
solar metallicity made our chemically consistent classical 
elliptical model redder at intermediate redshifts than the 
corresponding model for Zq. While we found good agree- 
ment between our chemically consistent galaxy models 
and the bulk of the HDF data over the redshift range 
from z > 4 through z > 0, a number of galaxies remained 
with colors bluer than those of our bluest undisturbed Sd 
model, that is described by a constant SFR and fairly 
low metallicity at intermediate and high redshifts, as well 
as quite a few galaxies with significantly redder colors in 
the redshift range 0.5 < z < 4 than is our reddest undis- 
turbed elliptical galaxy model, despite its high metallicity. 



This caused us to investigate to what extent starbursts 
can possibly explain the bluer-than-normal galaxies and 
to what extent post-starbursts might explain the redder- 
than-normal ones. 

In this paper, we studied the influence of starburst 
on the photometric evolution of galaxies at high red- 
shifts. Based on our evolutionary synthesis code GALEV, 
we calculated the spectrophotometric evolution of galax- 
ies bursting at redshifts 0.5, 1, 2, 3, and 4. To trans- 
late the time evolution of the models into a redshift 
evolution we used a standard cosmological model with 
Ho =70, rim = 0.3, r^A = 0.7 and an assumed redshift of 
galaxy formation Zf ~ 5. The star formation histories of 
undisturbed galaxy types E through Sd are defined by the 
requirement that models after a Hubble time agree with 
average observed properties of the respective galaxy types 
in terms of (i) broad band colors from UV through NIR, 
(ii) template UV and optical spectra, (iii) emission and 
absorption line strengths, (vi) gas content, and (v) metal 
abundances. Starbursts on top of these undisturbed mod- 
els are described by a sudden increase in the SFR followed 
by a decline on an e-folding timescale of 2.5 x 10* yr. Burst 
strengths are defined by the ratio of stellar mass formed 
out of the evolving available gas reservoir. In this defini- 
tion a burst of given strength produces more stars in a 
gas-rich galaxy or in a galaxy at a higher redshift than in 
a gas-poor galaxy or in a galaxy at a lower redshift. From 
the model spectra we calculated magnitudes and colors in 
the Johnson-Cousins and HST filter systems. Neither dust 
nor AGN contributions have been included in our models 
at the present stage. 

Following our models in their time evolution through 
starburst and post-starburst phases, we found the follow- 
ing results: 

— The blueing and brightening during a burst is a rela- 
tively short-term phenomenon lasting < 300—500 Myr. 
The reddening and dimming after the burst are much 
stronger, and the red and very red phases after the end 
of a burst can last very long, i.e. for many Gyr, unless 
a galaxy resumes SF from gas (re-)accreted after the 
end of the burst. 

— Bursts occurring at later phases, although transform- 
ing smaller absolute amounts of gas into stars, cause 
stronger blueing, since the galaxy itself is already rela- 
tively red, than bursts occurring at earlier stages, when 
the galaxy itself is still young and blue. On the other 
hand, the reddening after the burst is much stronger 
during the long evolution after early bursts than dur- 
ing the shorter time interval between a more recent 
burst and the present. 

Putting starburst and post-starburst evolution into a 
cosmological context and including the non-linear trans- 
formation between time and redshift, evolutionary effects, 
cosmological dimming and band shifts, as well as atten- 
uation by intergalactic neutral hydrogen, we observe the 
following: 
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— A burst of given strength starting at 
z — 4, 3, 2, 1, 0.5 makes a galaxy brighter in 
apparent me, nij and bluer than the undisturbed 
model until redshifts z - 2.8, 2.2, 1.5, 0.7, 0.3, 
respectively. 

— Bursts at high redshift, although producing larger 
amounts of stars, cause a very small blueing only and 
do so only in very small redshift intervals. 

— During their post-burst evolution, galaxies with bursts 
at high redshift get very red and remain so until low 
redshifts. 

We discussed our results in relation to the Extremely 
Red Objects (EROs) characterized by (R — K) > 5 and 
found in ample numbers at redshifts 0.5 ^ z < 2. We find 
that: 

— After bursts at redshifts between 2 and 4—5, our model 
galaxies, even without dust, naturally evolve into the 
color range of EROs in the course of their passive 
postburst evolution. The classical initial collapse and 
passive evolution elliptical model does not reach col- 
ors as red as those of the EROs, at least not for the 
progenitors of present-day average luminosity ellipti- 
cals. Particularly massive and luminous, and therefor 
particularly red, monolithic initial collapse ellipticals 
might reach ERO colors in the appropriate redshift 
range. They would then, however, be much more lu- 
minous at redshifts > 1 than any observed galaxy in 
the HDP. The K < 20 EROs require strong bursts on 
top of early-type spirals or the intrinsically brighter 
merger-induced starbursts as precursers. This fits well 
with the observation that ^ 40 % of all EROs show 
early-type morphologies, if we bear in mind that a ma- 
jor merger remnant will soon be transformed into an 
early-type galaxy in terms of morphology, as well as in 
terms of photometry, as our models show. 

The dust-free models presented here cannot, of course, 
describe the dusty-starburst fraction of ERO galaxies. 

A very first comparison of our models to HDF-N galax- 
ies with photometric redshifts from Sawicki et al. (1997) 
shows that almost all of the galaxies that were not ex- 
plained by our chemically consistent models for undis- 
turbed galaxies E . . . Sd are naturally explained by our 
starburst and post-starburst models: 

— HDP galaxies bluer than our bluest low-mctallicity, 
undisturbed Sd galaxy model, most of which are at low 
redshifts (z < 1) and a few also at higher redshifts, are 
explained well by ongoing starbursts. 

— HDP galaxies redder than our high-metallicity, undis- 
turbed, classical elliptical model, most of which are 
at redshifts around z ~^ 1, are well- matched by post- 
starburst models. The reddest galaxies around z ^ 1 
require models bursting at z > 2. 

Dust is not yet included in a self-consistent way in our 
models at the present stage. We have argued that the good 



agreement between the dust-free models and the data in- 
dicates that galaxies detected in UBVI are not strongly 
affected by dust. In dusty starbursts our models would, 
of course, severely underestimate the true burst strength. 
During post-starburst stages, observations of local galax- 
ies show that dust no longer plays an important role. Very 
shortly after a burst, galaxies apparently manage to get 
rid of most of their dust. 

The consistent inclusion of appropriate amounts of 
dust as well as a more extensive and quantitative in- 
terpretation of intermediate and high-redshift galaxies in 
terms of undisturbed, starburst, postburst galaxies, SPRs, 
mctallicity and mass evolution, etc., on the basis of our 
own photometric redshifts is in progress and will be the 
subject of a forthcoming paper. 
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